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We present new measurements of the inclusive forward- backward tt production asymmetry, Afb, 
and its dependence on several properties of the tt system. The measurements are performed with 
the full Tevatron data set recorded with the CDF II detector during pp collisions at y's = 1.96 TeV, 
corresponding to an integrated luminosity of 9.4 fb^ 1 . We measure the asymmetry using the rapidity 
difference Ay — y t — yi- Parton-level results are derived, yielding an inclusive asymmetry of 0.164± 
0.045 (stat + syst). We observe a linear dependence of Afb on the top-quark pair mass M t j 
and the rapidity difference |Aj/| at detector and parton levels. Assuming the standard model, 
the probabilities to observe the measured values or larger for the detector-level dependencies are 
7.4 x 10~ 3 and 14.7 x 10~ 3 for M ( j and |Aj/| respectively. Lastly, we study the dependence of the 
asymmetry on the transverse momentum of the ti system at the detector level. These results are 
consistent with previous lower-precision measurements and provide additional quantification of the 
functional dependencies of the asymmetry. 
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The creation of top quarks in pp collisions offers a 
unique test of pair-production in quantum chromody- 
namics (QCD) at very large momentum transfer as well 
as a promising potential avenue for the observation of 
new physical phenomena. Given the very large mass 
of the top quark, exotic processes may couple more 
strongly to top quarks than to the other known funda- 
mental particles, and possible hints of new interactions 
could be first observed in top-quark production. In par- 
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ticular, asymmetries in ti production could provide the 
first evidence of new interactions, such as ti production 
via a heavy axial color octet or a flavor-changing Z' 
boson, that might not be easily observed as excesses in 
the top quark production rate or as resonances in the 
ti invariant mass distribution. 

The CDF and DO collaborations have previously 
reported on forward-backward asymmetries (Ape) in 
pp — > ti production at y/s — 1.96 TeV at the Fer- 
milab Tevatron. In the standard model (SM), the ti 
production process is symmetric in production angle 
except for a 0(7%) charge asymmetry arising at next- 
to-leading order (NLO) and beyond in QCD pQ. Us- 
ing a sample corresponding to 5.3 fb _1 of integrated 
luminosity, CDF measured a parton-level asymmetry 
A FB = 0.158 ± 0.074 in the lepton+jets decay chan- 
nel (ti (W+b)(W-b) -> (l+v)(qq')bb 0), and very 
good agreement was found by the DO measurement 
Afb = 0.196 ± 0.065 0] in a lepton+jets sample cor- 
responding to 5.4 fb _ . CDF and DO have also per- 
formed simple differential measurements using two bins 
each in the top-antitop rapidity difference \Ay\ and the 
top-antitop invariant mass M t j. The two experiments 
agreed on a large |Ay| dependence. CDF also saw a 
large M t i dependence, and while that observed at DO 
was smaller, the CDF and DO results were statistically 
consistent. One of the aims of this paper is to clarify 
the | Ay | and M t i dependence of the asymmetry using 
the full CDF data set. 

The 5 fb _1 results have stimulated new theoretical 
work, both within and outside the context of the SM. 
The SM calculation has been improved by calculations 
of electroweak processes that contribute to the asymme- 
try, studies of the choice of renormalization scale, and 
progress on a next-to-next-to-leading order (NNLO) 
calculation of the asymmetry [5--9J . The new calcu- 
lations result in a small increase in the expected asym- 
metry, but not enough to resolve the tension with ob- 
servation. Other work has focused on the dependence 
of the asymmetry on the transverse momentum of the 
ti system |10j . on which we report here. 

A number of speculative papers invoke new interac- 
tions in the top sector [TT] to explain the large asym- 
metry. In one class of models, ti pairs can be pro- 
duced via new axial s-channel particles arising from 
extended gauge symmetries or extra dimensions. For 
these models, the asymmetry is caused by interference 
between the new s-channel mediator and the SM gluon. 
In other models, light t-channel particles with flavor- 
violating couplings create an asymmetry via a u, d — > t 
flavor change into the forward Rutherford-scattering 
peak. All potential models of new interactions must ac- 
commodate the apparent consistency of the measured 
cross section and M t i spectrum with the SM predic- 
tions. Tevatron and LHC searches for related phenom- 
ena, such as di-jet resonances, same-sign tops, and other 
exotic processes, can provide additional experimental 



limits on potential models. Measurements by the LHC 
experiments of the top-quark charge asymmetry Aq 1 
an observable that is distinct from Apb but correlated 
with it, have found no significant disagreement with the 
SM [12] : however, any observable effect at the LHC is 
expected to be small, and the nature of the relationship 
between Afb and Ac is model-dependent [13] . A more 
precise measurement of the Tevatron forward-backward 
asymmetry and its mass and rapidity dependence may 
help untangle the potential new physics sources for A F b 
from the standard model and from each other. 

This paper reports on a study of the asymmetry in the 
lepton+jets topology, with several new features com- 
pared to the previous CDF analysis in this channel [2J. 
We use the complete Tevatron Run II data set with an 
integrated luminosity of 9.4 fb~ . We additionally ex- 
pand the online event selection by including events trig- 
gered by large missing transverse energy and multiple 
hadronic jets, increasing the total data set by approx- 
imately 30% beyond what is gained by the increase in 
luminosity In total, the number of candidate events in 
this analysis is more than twice the number of events 
used in Ref. [5] . An improved NLO Monte Carlo gener- 
ator is used to describe the predicted ti signal, and we 
also add small corrections reflecting new results on the 
electroweak contributions to the asymmetry j5H2]- Fi- 
nally, parton-level shape corrections utilize an improved 
algorithm which yields multi-binned parton-level mea- 
surements of the rapidity and mass dependence of the 
asymmetry. We also study the dependence of the asym- 
metry on the ti transverse momentum, p", showing that 
the modeling of this quanity is robust, and that the ex- 
cess asymmetry above the SM prediction is consistent 
with being independent of p". 



II. EXPECTED ASYMMETRIES AND MONTE 
CARLO MODELS 

The asymmetry is measured using the frame- 
invariant difference of the t and i rapidities, Ay = 
Vt — Vti where the rapidity y is given by 

with E being the total top-quark energy and p z being 
the component of the top-quark momentum along the 
beam axis as measured in the detector rest frame. Ay 
is invariant to boosts along the beamline, and in the 
limit where the transverse momentum of the ti system 
is small, the forward-backward asymmetry 

N(Ay > 0) - N(Ay < 0) 
FB N(Ay > 0) + N(Ay < 0) [ ' 
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is identical to the asymmetry in the top-quark pro- 
duction angle in the tt rest frame. The standard 
model predictions for the top-quark asymmetry refer- 
enced in this paper are based on the NLO event gener- 
ator powheg [14], using the CTEQ6.1M set of parton- 
distribution functions (PDFs), validated by comparing 
POWHEG to the NLO generator MC@NLO [15] as well 
as the NLO calculation of mcfm [TBI . We find good 
consistency overall, as shown in Table ffl [17] . Sources of 
asymmetry from electroweak processes in the standard 
model that are not included in the Monte Carlo calcula- 
tions [MZ] lead to an overall increase of the asymmetry 
by a factor of 26% of the QCD expectation, which is 
included in all the predictions shown in Table [T] All pre- 
dicted asymmetries and Ay distributions in this paper 
contain this multiplicative 26% electroweak correction 
to POWHEG. Because the electroweak tt production pro- 
cesses that contribute to the asymmetry correspond to 
similar QCD processes (after replacing the gluon with 
a photon or Z boson), the electroweak asymmetry is 
assumed to have a dependence on M t j and Ay similar 
to the QCD asymmetry. Thus, the correction increases 
the QCD predictions by a constant factor of 26% for 
all values of M t i and Ay. Following Ref. [18] . we in- 
clude a 30% uncertainty on all theoretical predictions 
for the SM asymmetry due to the choice of renormal- 
ization scale. 

To test the analysis methodology in the case of a large 
asymmetry, we study two models in which an asymme- 
try is generated by the interference of the gluon with 
massive axial color-octet particles. Each provides a 
reasonable approximation of the observed data in pre- 
senting a large, positive forward-backward asymmetry, 
while also being comparable to the data in other impor- 
tant variables such as the tt invariant mass, M t j. 

The first model, Octet A, contains an axigluon with 
a mass of 2 TeV/c 2 . This hypothetical particle is mas- 
sive enough that the pole is observed as only a small 
excess in the tail of the M t i spectrum, but it creates 
an asymmetry via the interference between the off-shell 
axigluon and the SM gluon. The couplings are tuned 
(9v(q) = 9v(t) = 0, g A (q) = 3, g A (t) = -3, where q 
refers to light-quark couplings and t to top-quark cou- 
plings) to produce a parton-level asymmetry consistent 
with the measurement in Ref. [2J. The second model, 
Octet B, contains an axigluon with the same couplings, 
but a smaller mass of 1.8 TeV/c 2 . This model produces 
a larger excess in the tail of the M t j spectrum and an 
even larger asymmetry than Octet A, allowing the mea- 
surement procedure to be tested in a regime with a very 
large asymmetry. 

Both models are simulated using the leading order 
(LO) MADGRAPH [HI [20] Monte Carlo generator and 
are hadronized with PYTHIA [55] before being passed to 
the CDF detector simulation and reconstruction soft- 
ware. We emphasize that these are not hypotheses - the 
physical applicability of these models is, in fact, quite 



constrained by tt resonance searches at the LHC |21) . 
Rather, these models are use as controlled inputs to 
study the performance of the analysis in the presence 
of large asymmetries. Further information about these 
models can be found in Ref. [2]. 



III. MEASUREMENT STRATEGY 

The analysis takes place in several steps. We first 
consider the asymmetry observed at the reconstruction 
level in all selected events. Next, to study the asym- 
metry for a pure sample of tt events as recorded in 
the detector, the calculated non-tt background contri- 
bution is subtracted and the appropriate systematic un- 
certainties related to the background prediction are ap- 
plied. Finally, to study the asymmetry at the parton 
level, corrections are applied for the event reconstruc- 
tion and detector acceptance, along with appropriate 
systematic uncertainties on the signal modeling. The 
reconstruction- and background-subtracted-level mea- 
surements have the advantage of fewer assumptions, 
while the parton-level measurement allows direct com- 
parison to theory predictions. 

After reviewing the event selection and reconstruc- 
tion in Sec. |IV| we describe the various steps of the 
correction procedure in detail and apply them to the 
Ay distribution and the inclusive Afb measurement in 
Sec. |V] In Sec. [V| and Sec. [VIlJ we study the depen- 
dence of the asymmetry on \Ay\ and M t j, Aps(\Ay\) 
and ApB{M t t) respectively, at all three stages of cor- 



rection, and Sec. |VIII| discusses the significance of dis- 
crepancies observed in these dependencies between the 
data and the SM. Section [IX] discusses the dependence 
of the asymmetry on the tt transverse momentum. 



IV. DETECTOR, EVENT SELECTION, AND 
RECONSTRUCTION 

The data sample corresponds to an integrated lumi- 
nosity of 9.4 fb _1 recorded with the CDF II detector 
during pp collisions at 1.96 TeV. CDF II is a general 
purpose, azimuthally and forward-backward symmetric 
magnetic spectrometer with calorimeters and muon de- 
tectors [53]. Charged particle trajectories are measured 
with a silicon-microstrip detector surrounded by a large 
open-cell drift chamber, both within a 1.4 T solenoidal 
magnetic field. The solenoid is surrounded by pointing- 
tower-geometry electromagnetic and hadronic calorime- 
ters for the measurement of particle energies and miss- 
ing energy reconstruction. Surrounding the calorime- 
ters, scintillators and proportional chambers provide 
muon identification. We use a cylindrical coordinate 
system with the origin at the center of the detector and 
the z-axis along the direction of the proton beam [24] . 
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TABLE I: Parton-level asymmetry predictions of powheg, mc@nlo, and mcfm after applying electroweak corrections. 







MC@NLO 


POWHEG 


MCFM 




Inclusive 


0.067 ±0.020 


0.066 ±0.020 


0.073 ±0.022 




|Ay| < 1 
|Ay| > 1 


0.047 ±0.014 
0.130 ±0.039 


0.043 ±0.013 
0.139 ± 0.042 


0.049 ±0.015 
0.150 ±0.045 


M a 


< 450 GeV/c 2 
> 450 GeV/c 2 


0.054 ±0.016 
0.089 ±0.027 


0.047 ± 0.014 
0.100 ±0.030 


0.050 ±0.015 
0.110 ±0.033 



This measurement selects tt candidate events in 
the lepton±jets topology, where one top quark decays 
semileptonically (t — > Wb — > Ivb) and the other hadron- 
ically it —> Wb —> qq'b). We detect the lepton and 
hadronization-induced jets. The presence of missing 
transverse energy (Et) [21] is used to infer the pas- 
sage of a neutrino through the detector. Detector read- 
out is initiated in one of two ways: either by indica- 
tions of a high-momentum lepton (electron or muon) 
in the central portion of the detector or by events 
with indications of large Et and at least two ener- 
getic jets. Events collected in the second manner, in 
which we require the presence of muon candidates re- 
constructed offline, make up the "loose muon" sample, 
a new addition compared to the previous version of this 
analysis. After offline event reconstruction, we require 
that all candidate events contain exactly one electron 
or muon with Et(pt) > 20 GeV(GeV/c), as well as 
four or more hadronic jets with Ex > 20 GeV and 
\r)\ < 2.0. Jets are reconstructed using a cone algo- 
rithm with AR = \J A<j) 2 ± A?/ 2 < 0.4, and calorime- 
ter signals are corrected for various detector and mea- 
surement effects as described in Ref. [23]. We require 
Et > 20 GeV, consistent with the presence of an unde- 
tected neutrino. We finally require that Ht, the scalar 
sum of the transverse energy of the lepton, jets, and Et, 
be Ht > 220 GeV. This requirement reduces the back- 
grounds by 17% while accepting 97% of signal events. 
The SECVTX algorithm [25] is used to identify b jets 
by searching for displaced decay vertices within the jet 
cones, and at least one jet in each event must contain 
such a "6 tag". The coverage of the tracking detector 
limits the acceptance for jets with identified b tags to 

M<i. 

The sample passing this selection, including the 6-tag 
requirement, contains 2653 candidate events. The esti- 
mated non-tt background in the data sample is 530±124 
events. The predominant background source is QCD- 
induced W^±multi-parton events containing either b- 
tagged heavy-flavor jets or erroneously tagged light- 
flavor jets. These events are modeled with the ALPGEN 
Monte Carlo generator [27] , with the normalizations de- 
termined by tagging efficiencies, mis-tagging rates, and 
other measurements in the data. QCD multi-jet ("Non- 
W") events containing mis-measured Et and jets that 



are mis-identified as leptons are modeled using real data 
events with lepton candidates that are rejected by the 
lepton identification requirements. This background, 
which is the most difficult to model properly, is also 
the one that is most efficiently suppressed by the Ht 
requirement, which reduces it by approximately 30%. 
Small backgrounds from electroweak processes (WW, 
WZ, single-top) are estimated using Monte Carlo gen- 
erators. The expected background contributions from 
each source are given in Table |H] We note that there are 
correlations among the various sources of uncertainty 
for the different background components, so that the 
total background uncertainty is not a simple sum in 
quadrature of the uncertainties on the individual back- 
ground normalizations. Further information about the 
background modeling and event selection can be found 
in Ref. [25] . 

TABLE II: Expected contributions of the various back- 
ground sources to the selected data. 



Background source 


Number of events 


W+UF 


256 ± 83 


W+LF 


102 ± 32 


Non-VK 


97 ± 50 


Single top 


35 ± 3 


Diboson 


21 ± 3 


Z±jets 


19 ± 3 


Total background 


530 ± 124 


ti (7.4 pb) 


2186 ± 314 


Total prediction 


2716 ± 339 


Data 


2653 



The reconstruction of the tt kinematics employs the 
measured momenta of the lepton and the four leading 

jets in the event, along with the measured Et- The cal- 
culation of the tt four- vectors uses a x 2 -based fit of the 
lepton and jet kinematic properties to the tt hypothesis. 
Each of the possible jet-to-parton assignments is evalu- 
ated according to its consistency with resulting from the 
decay of a pair of top quarks. Two of the observed jets 
are required to be consistent with being decay products 
of a W boson, while the lepton and Et must be con- 
sistent with another W boson. Each W boson, when 
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paired with one of the remaining (6) jets, is checked for 
consistency with having resulted from a top-quark de- 
cay. The lepton momentum, ^r, and jet energies are 
allowed to float within their experimental uncertainties, 
and we apply the constraints that Mw = 80.4 GeV/c 2 , 
M t = 172.5 GeV/c 2 , and any ^-tagged jets must be as- 
sociated with b partons. The jet-to-parton assignment 
that best matches these requirements is chosen to define 
the parent top quarks in each event. 

This algorithm has been studied and validated in 
many precision top-quark-property analyses, including 
mass measurements |29j . which remove the top-quark 
mass constraint, and property measurements that do 
make use of the mass constraint |30j . The top- and 
antitop-quark four-vectors determined from this proce- 
dure are used to find the rapidities of the quarks and the 
Ay = yt — Hi variable used for the asymmetry analysis, 
with the charges of the reconstructed top quarks being 
fixed by the observed lepton charge. In the Appendix, 
we discuss a high-precision test of the lepton-charge de- 
termination in a large control sample with the goal of 
verifying that the lepton charge is correctly assigned. 

The validity of the analysis is checked at all stages by 
comparison to a standard model prediction created us- 
ing the powheg tt model, the lepton+jets background 
model described above, and a full simulation of the 
CDF II detector. Figure [T] shows the rapidity distribu- 
tion for the hadronically-decaying top or antitop quark. 
In the measurement of the asymmetry, the observed 
lepton charge is used to determine whether each entry 
in this distribution corresponds to a top quark or an 
antitop quark, and this rapidity is combined with the 
rapidity of the leptonically decaying quark to calculate 
Ay for each event. In Fig. [T] and all that follow, the tt 
signal prediction is scaled such that the total signal nor- 
malization, when added to the background prediction in 
Table ITT1 totals number of observed events. 
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FIG. 1: The rapidity of the hadronically-decaying top or 
antitop quark. 



Figure [2] shows a comparison of the data to the pre- 



diction for the invariant mass of the tt system, M t f, 
there is good agreement. In the previous CDF analy- 
sis [2], the forward-backward asymmetry was found to 
have a large dependence on this variable. In Sec 
we report a new measurement of this dependence. 
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FIG. 2: Reconstructed invariant mass of the tt system. The 
last bin contains overflow events. 



The transverse momentum of the tt system, p", pro- 
vides a sensitive test of the reconstruction and model- 
ing, particularly at low momenta, where both the pre- 
diction and the reconstruction are challenged by the ad- 
dition of soft gluon radiation external to the tt system. 
In the reconstructed data, we find good agreement with 
the sum of the background prediction and the NLO tt 
model, as shown in Fig. [3] The tt forward-backward 
asymmetry can have a significant dependence, and 
we discuss the expected and measured asymmetry as a 
function of this variable in Sec. IIXI 
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FIG. 3: Reconstructed of the tt system. The last bin 
contains overflow events. 



We also consider a wide range of other variables, a 
selection of which are shown here, to validate the re- 
construction algorithm and the modeling of the data 
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FIG. 4: (a) The number of observed jets and |(b)| the number of jets with b tags in the data compared to the signal plus 



background model. The last bin contains overflow events. 



set. In Fig. [4] we show the distributions of the number 
of jets and number of b tags in events passing the se- 
lection requirements. Figure |4(b)| also includes events 
containing no fe-tagged jets, which are not part of the 
final sample of candidate events but provide an impor- 
tant check on the modeling of the 6-tagging algorithm. 
Figure [5] shows the transverse energy of the most ener- 
getic jet and the transverse momentum of the lepton, 
while Fig. [6] shows the distribution of the reconstructed 
$t and Ht- All distributions exhibit good agreement 
between the observed data and the model expectations. 



V. THE INCLUSIVE ASYMMETRY 
A. Ay in the reconstructed data 



effect. 



TABLE III: Measured reconstruction-level asymmetries in 
Ay compared to the values measured in the previous CDF 
analysis fJJ, as well as the predicted asymmetries for the 
signal and background contributions. 



Predicted Afb 


SM tt 


0.033 ±0.011 


Backgrounds 


-0.034 ±0.013 


Total prediction 


0.020 ±0.012 




Observed j4fb± stat 




9.4 fb" 1 5.3 fb" 1 


All data 


0.063 ±0.019 0.057 ±0.028 


Positive leptons 


0.072 ±0.028 0.067 ±0.040 


Negative leptons 


0.055 ± 0.027 0.048 ± 0.039 



We first consider the reconstructed Ay distribution 
and its asymmetry as defined in Eq. The Ay dis- 
tribution is shown in Fig. [7J compared to prediction 
for the background plus the POWHEG tt model. Those 
bins with Ay > contain data points that are consis- 
tently higher than the prediction, while in the bins with 
Ay < 0, the data is consistently below the prediction. 
This results in an inclusive reconstructed asymmetry 
of ApQ — 0.063 ± 0.019, compared to a prediction of 
0.020 ± 0.012. The uncertainty on the data measure- 



ment is statistical only. Table III summarizes the recon- 
structed asymmetry values, with events split according 
to the charge of the identified lepton, and also reports 
the results of Ref. [5] for comparison. The uncertain- 
ties scale as expected from the previous analysis accord- 
ing to the increase in the number of candidate events. 
When the sample is separated according to the charge 
of the lepton, the asymmetries are equal within uncer- 
tainties, as would be expected from a CP-conserving 



B. Subtracting the background contributions 

Approximately 20% of the selected data set is com- 
posed of events originating from various background 
sources. We remove the effect of these events by sub- 
tracting the predicted background contribution from 
each bin of the reconstructed distribution. This 
background-subtraction procedure introduces a system- 
atic uncertainty, which is added in quadrature to the 
statistical uncertainty for all background-subtracted re- 
sults in this paper. 

The total systematic uncertainty on quantities deter- 
mined from the background-subtracted data arises from 
an uncertainty on the background normalization and an 
uncertainty on the background shape, and is given by 
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FIG. 6: 

data compared to the signal plus background model. The last bin contains overflow events. 



(3) 

To derive this uncertainty, we start with a total pre- 
diction containing n components (one signal and n — 1 
background sources), with each component i having an 
asymmetry and contributing A/j events. This leads 
to a total asymmetry for the prediction of 



E AiNi £ AiNi 
Atot = ^4 = ■ (4) 

E Ni Ntot 

i=l 



For the iih component, we let OAi and be the un- 
certainties on the asymmetry and the normalization re- 
spectively. These uncertainties can be propagated in 
the usual way by calculating derivatives and adding in 
quadrature, leading to the term within the summation 
in Eq. ([3| . For the uncertainty due to background sub- 
traction, the summation runs over the n— 1 background 
components. We also include an overall uncertainty 
related to the assumptions associated with the back- 
ground modeling, resulting in the last term in Eq. 

For <TN i , we use the uncertainty on the predicted nor- 
malization of each background component, as listed in 
Table [TTJ The uncertainty due to the finite size of the 
Monte Carlo sample used to model a given background 
component is included as a At . For the inclusive Afb 
measurement, the summation term in Eq. (13]) results 
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FIG. 7: (top) The reconstructed Ay distribution and the 
inclusive reconstruction-level asymmetry, compared to the 
prediction of the signal and background model, (bottom) 
The difference between the data and prediction divided by 
the prediction. 



FIG. 8: (top) The observed background-subtracted Ay dis- 
tribution compared to the SM prediction. Error bars include 
both statistical and background-related systematic uncer- 
tainties, (bottom) The difference between the data and pre- 
diction divided by the prediction. 



in a total uncertainty of 0.008. The systematic uncer- 
tainty (Tyi bkg on the overall background shape is esti- 
mated by imposing reasonable variations on the back- 
ground model and determining the effect on the mea- 
sured asymmetry, with the last term in Eq. ([3| yielding 
an uncertainty of 0.002 for the inclusive Ape result. In 
total, the sum of the systematic contributions to the 
uncertainty is small compared to the statistical uncer- 
tainty. 

The Ay distribution after background subtraction is 
shown in Fig. [8j Although the total background pre- 
diction has only a small intrinsic asymmetry, the re- 
moval of the backgrounds decreases the dilution of the ti 
signal asymmetry, resulting in an observed asymmetry 
in the background-subtracted sample of 0.087 ± 0.026 
(stat+syst), compared to the powheg prediction of 
0.033 ±0.011. 



C. Correction to the parton level 

The background-subtracted results provide a mea- 
surement of the asymmetry in a sample composed 
purely of ti events. However, these results are not 
directly comparable to theoretical predictions because 
they include the effects of the limited acceptance and 
resolution of the detector. We correct for these effects 
so as to provide parton-level results that can be directly 
compared to theory predictions. 

If the true parton-level binned distribution of a 
particular variable is given by n part on, then, after 
background subtraction, we will observe Hbkg.sub. = 
SAn par ton, where the diagonal matrix A encodes the 
effect of the detector acceptance and selection require- 
ments, while the response matrix S describes the bin- 



to-bin migration that occurs in events passing the se- 
lection due to the limited resolution of the detector and 
ti reconstruction algorithm. To recover the parton-level 
distribution, the effects of S and A must be reversed. 

The 5.3 fb _1 CDF analysis [5] used simple matrix in- 
version ("unfolding") to perform the correction to the 
parton level. While effective, this technique was limited 
in its application because matrix inversion tends to en- 
hance statistical fluctuations, which makes it reliable 
only in densely populated distributions. This limited 
the previous analysis to the extent that the determi- 
nation the functional dependencies of the asymmetry 
could only use two bins of |Ay| and M t i. In this paper, 
we employ a new algorithm, also based on matrix inver- 
sion but more sophisticated in application, to measure 
more finely-binned parton-level distributions, resulting 
in a more robust measurement of the functional depen- 
dence of A FB on | Ay | and M t j at the parton level. 

We first consider S, correcting for the finite resolution 
of the detector using a regularized unfolding algorithm 
based on Singular Value Decomposition (SVD) [3"Tll3"2"] . 
We model the bin-to-bin migration caused by the de- 
tector and reconstruction using powheg. The matrix 
S in Ay from powheg is represented graphically in 
Fig. [9] Along each row, the box area is proportional 
to the probability that each possible measured value 
Ay mcas is observed in events with a given true rapid- 
ity difference Ay truc . The matrix population clusters 
along the diagonal where Ay mcas = Ay true and is ap- 
proximately symmetric, showing no large biases in the 
Ay reconstruction. Before inverting the matrix S and 
applying it to the background-subtracted data, a regu- 
larization term is introduced to prevent statistical fluc- 
tuations from dominating the correction procedure. It 
is this smoothing via regularization that allows an in- 
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FIG. 9: Detector response in Ay as modeled by powheg, 
showing the true value of Ay as a function of the measured 
value for all events passing the selection criteria. The size 
of each rectangle is proportional to the number of entries in 
that bin. 



crease in the number of bins in the parton-level dis- 
tributions compared to the previous analysis. Details 
regarding how the regularization term is included are 
given in Ref. [31], but in essence, a term y/rC, where 
C is the second-derivative matrix, 



C = 



o 
1-2 1 
1-21 



V 



1 -2 
1 



(5) 



-1/ 



is added to the matrix equation relating n pa rton to 
^bkg.sub. so as to impose the a priori condition that 
the parton-level solution should be smooth. The min- 
imization of this regularization term (which is, in ef- 
fect, a minimization of the second derivatives of the 
parton-level result) is used as a constraint when the 
background-subtracted data nbkg.sub. is multiplied by 
the inverse of the response matrix S. The value of r de- 
fines how strongly the regularization condition affects 
the result and is determined using the methods recom- 
mended in Ref. [31] • If the smoothing is too strong, 
large biases can result, while if it is too weak, statisti- 
cal precision is lost; to guard against these possibilities, 
simulated experiments are performed as described be- 
low using various values for r in order to verify that the 
chosen strength yields both good statistical precision 
and minimal residual bias. 

In the second step of the parton-level correction pro- 
cedure, we account for events that are unobserved due 
to limited acceptance. The acceptance in each bin is de- 




FIG. 10: Acceptance as a function of Ay as modeled by 

POWHEG. 



these acceptances are applied to the data as an inverse- 
multiplicative correction to each bin. The acceptance is 
asymmetric in Ay, with backwards events passing the 
selection requirements more often than forward events. 
This effect is related to the dep endence of the asym- 

Large in a given 



IX 



metry that is discussed in Sec. 
event leads to tt decay products that also have large 
Pt, and thus events with large pass the selection re- 
quirements m ore o ften than events with small . As is 
shown in Sec. IX high-p" events are also predicted by 



powheg (and various other SM calculations) to have 
a negative asymmetry. The result is that events with 
a negative asymmetry are more likely to fulfill the se- 
lection requirements, leading to the asymmetric accep- 



tance distribution in Fig. 10 



rived from the powheg model, as shown in Fig. 10 and 



The SVD unsmearing and bin-by-bin acceptance cor- 
rection have similarly-sized impact on the final result. 
Both of the corrections lead to an increase in the asym- 
metry: the unsmearing increases App, by removing dilu- 
tion from mis-reconstructed events, which tend to have 
no asymmetry, and the acceptance correction increases 
Afb because of the asymmetric acceptance shown in 
Fig. [101 

The combination of these two parts of the correction 
procedure allows the determination of the parton-level 
distribution of Ay, which is reported as a differential 
cross section. This algorithm is tested in various sim- 
ulated tt samples, including standard model powheg 
and the non-SM samples Octet A and Octet B. Analyz- 
ing these samples as if they are data, we measure the 
bias in the comparison of derived parton-level results to 
the true values in the generated samples. The powheg 
results are self-consistent to better than 1%, and, be- 
cause the NLO standard model is assumed a priori to 
be the correct description of the underlying physics and 
is used to model the acceptance and detector response, 
any biases observed in this case are included as system- 
atic uncertainties, as described below. 
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TABLE IV: Average parton-level asymmetry values in 10 000 simulated experiments with Octet A. 



|Ay| 
Inclusive 
0.0 < |Ay| < 0.5 
0.5 < jAyj < 1.0 
1.0 < \Ay\ < 1.5 
\Ay\ > 1.5 



Average measured Afb 
0.162 
0.056 
0.180 
0.316 
0.434 



Average uncertainty 
0.039 
0.035 
0.055 
0.078 
0.128 



True Afb 
0T56 
0.052 
0.158 
0.295 
0.468 
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FIG. 11: Results from simulated parton-level Ay measure- 
ments based on Octet A. The data points show the central 
values for the simulated results, with the error bars repre- 
senting the la spread of the results. 



TABLE V: Systematic uncertainties on the parton level Afb 
measurement. 



Source 


Uncertainty 


Background shape 


0.014 


Background normalization 


0.013 


Hadronization 


0.010 


Jet energy scale 


0.005 


Initial- and final-state radiation 


0.005 


Correction procedure 


0.003 


Color reconncction 


0.001 


Parton-distribution functions 


0.001 


Total systematic uncertainty 


0.023 


Statistical uncertainty 


0.039 


Total uncertainty 


0.045 



In the octet models, the derived distributions track 
the generator truth predictions well, but small biases 
(generally less than 15% of the measured asymmetries) 
are observed in some of the differential asymmetry val- 
ues. An example of the average corrected distribution 
across a set of 10 000 simulated experiments is shown in 
Fig. 



11 for Octet A, with the asymmetry as a function 



of | Ay | for these simulated experiments summarized in 
Table |IV) We do not attempt to correct the biases seen 
in the non-SM models or include them in the uncer- 
tainty because there is no reason to believe that these 
specific octet models actually represent the real under- 
lying physics - these models exhibit small but significant 
discrepancies with the data in the M t j spectrum, a vari- 
able that has a significant effect on the tt reconstruction, 
and thus the detector response matrix. In light of this 
model-dependence, we emphasize that the parton-level 
results need to be interpreted with some caution in re- 
lation to models that differ significantly from the NLO 
standard model. 

Because the resolution corrections can cause migra- 
tion of events across bins, the populations in the final 
parton-level distributions are correlated. In all binned 
parton-level distributions, the error bars on a given bin 
correspond to the uncertainty in the contents of that 
bin, but they are not independent of the uncertainties 
corresponding to other bins in the distribution. When 
we calculate derived quantities such as j4fb, we use the 
covariance matrix associated with the unsmearing pro- 
cedure to propagate the uncertainties correctly. 

Several sources of systematic uncertainty must be ac- 
counted for when applying the correction procedure. 
In addition to uncertainties on the size and shape of 
the background prediction, there are also uncertain- 
ties related to the signal Monte Carlo sample used to 
model the acceptance and detector response, including 
the size of the jet energy scale corrections, the amount 
of initial- and final-state radiation, the models used for 
color reconnection and hadronization, and the under- 
lying parton-distribution functions. We evaluate these 
uncertainties by repeating the measurement after mak- 
ing reasonable variations to the assumptions that are 
used when modeling the detector response. We also 
include a systematic uncertainty for the correction al- 
gorithm itself, taking the difference between the true 
value in powheg and the average result from the sim- 
ulated experiments based on POWHEG described above 
as the uncertainty resulting from the correction pro- 
cedure. The systematic uncertainties on the inclusive 
Afb measurement are shown in Table [Vj and the total 
systematic uncertainty is found to be small compared 
to the statistical uncertainty. 
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FIG. 12: (top) The differential cross section da/d(Ay) as 
measured in the data after correction to the parton level 
compared to the SM prediction. Uncertainties include both 
statistical and systematic contributions and are correlated 
between bins, (bottom) The difference between the data 
and prediction divided by the prediction. 



Applying the correction procedure to the data of 
Fig. [8] yields the distribution shown in Fig. [l2j where 
the measured result is compared to the SM powheg 
prediction. Both the prediction and the observed data 
distributions are scaled to a total cross section of 7.4 pb, 
so that Fig. 12 shows the differential cross section for 

We measure an in- 



tt production as a function of Ay. 
elusive parton-level asymmetry of 0.164 ± 0.039(stat) ± 
0.023(syst) = 0.164 ± 0.045. At the parton level, the 
observed inclusive asymmetry is non-zero with a sig- 
nificance of 3.6a and exceeds the NLO prediction of 
powheg by 2.0ct, where we have included a 30% uncer- 
tainty on the prediction. 



VI. THE DEPENDENCE OF THE 
ASYMMETRY ON \Ay\ 

The dependence of A FB on the rapidity difference 
| Ay | was studied in the 5 fb -1 analyses OH], but with 
only two bins of |Ay|. The CDF and DO results were 
consistent and showed a rise of ^4fb with increasing 
\Ay\. We perform a more detailed study of the rapidity 
dependence of A F b using the full data set and improved 
analysis techniques. 

The forward-backward asymmetry as a function of 
| Ay | at the reconstruction level can be derived from the 
data shown in Fig. [7] according to 



AwB(\Ay\) 



N F (\Ay\) - N B (\Ay\) 
N F (\Ay\)+N B (\Ay\y 



(6) 



where NpQ Ay\) is the number of events in a given | Ay| 
bin with Ay > and A^b(|Aj/|) is the number of events 
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FIG. 13: The reconstruction-level forward-backward asym- 
metry as a function of | Ay\ with a best-fit line superimposed. 
The errors on the data are statistical, and the shaded region 
represents the uncertainty on the slope of the prediction. 



in the corresponding \Ay\ bin with Ay < 0. Figure 13 
shows Afb(|A2/|) in four bins of \Ay\, with the mea- 
sured values and their uncertainties listed in Table IVII 
To quantify the behavior in a simple way, we assume 
a linear relationship, which provides a good approxi- 
mation of both the data and the powheg prediction 
(see also Ref . [33J ) ■ Additionally, the effect of interest - 
the dependence of A F b on \Ay\ - can be described en- 
tirely by the slope of the resulting line. The intercept 
yields no information about the dependence j4fb(|Aj/|), 
and in fact, we expect the intercept to be consistent 
with zero within the statistical sensitivity. Any theory 
that predicts a continuous and differentiable Ay dis- 
tribution must have Afb(|Aj/| = 0) = regardless of 
the size of the inclusive asymmetry. The slope a& y of 
the line does not correspond to a specific parameter of 
any particular theory, but provides a quantitative com- 
parison of the | Ay | behavior of the asymmetry in the 
data and prediction. The measurements of Afb(|Aj/|) 
in data at the reconstruction level are well-fit by a line 
with a x 2 P er degree of freedom of 1.1/2 and a slope 
a A y — (14.3 ± 4.3) x 10~ 2 , a rapidity dependence that 
is non-zero with significance in excess of 3c. The pre- 
dicted slope from powheg and the background model 
is (5.6 ± 1.4) x 10~ 2 . 

The behavior of the asymmetry as a function of 
| Ay | is also measured after the removal of the back- 
ground contribution as described previously. Figure 14 
shows the distribution v4pB(|Ay|) for the background- 
subtracted ti signal, with the measured values sum- 
marized in Table VII Systematic uncertainties on the 



background-subtraction procedure are included in the 
error bars. The data measurements and the predictions 
are well-fitted by the linear assumption, with an ob- 
served slope of aA y = (18.2 ± 5.7) x 10~ 2 that exceeds 
the prediction of (6.6±1.6) x 10 -2 by approximately 2a. 
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TABLE VI: The asymmetry at the reconstructed level as measured in the data, compared to the SM tt plus background 
expectation, as a function of [Ay|. 





Data 


SM tt + Bkg. 


\Ay\ 


^4fb ± stat 




0.0 - 0.5 


0.016 ± 0.028 


0.001 ± 0.005 


0.5 - 1.0 


0.055 ± 0.035 


0.020 ± 0.012 


1.0 - 1.5 


0.186 ± 0.049 


0.050 ± 0.021 


> 1.5 


0.206 ± 0.085 


0.109 ± 0.030 



TABLE VII: The asymmetry at the background-subtracted level as measured in the data, compared to the SM tt expectation, 
as a function of jAyj. 





Data 


SM tt 


|Aj/| 


Afb ± (stat+syst) 


A FB 


0.0 - 0.5 


0.027 ± 0.034 


0.009 ± 0.005 


0.5 - 1.0 


0.086 ± 0.045 


0.040 ± 0.014 


1.0 - 1.5 


0.246 ± 0.063 


0.074 ± 0.026 


> 1.5 


0.254 ± 0.124 


0.113 ± 0.039 



DQ 

< 




FIG. 14: The background-subtracted asymmetry as a func- 
tion of |Ay| with a best-fit line superimposed. Error bars 
include both statistical and background-related systematic 
uncertainties. The shaded region represents the theoretical 
uncertainty on the slope of the prediction. 




Pa rton- Level |Ay| 



FIG. 15: The parton-level forward- backward asymmetry as 
a function of \Ay\ with a best-fit line superimposed. Un- 
certainties are correlated and include both statistical and 
systematic contributions. The shaded region represents the 
theoretical uncertainty on the slope of the prediction. 



The observed slope is larger than at the reconstruction 
level owing to the removal of the background dilution, 
with the significance of the difference relative to the 
standard model staying approximately the same. 

The | Ay | dependence of the asymmetry at the par- 
ton level can be derived from Fig. [12] by comparing the 
forward and backward bins corresponding to a given 
value of \Ay\. This parton-level A FB (|Ay|) distribution 
is shown in Fig. |15| wit h the asymmetries in each bin 
also listed in Table VIII A linear fit to the parton-level 



results yields a slope a Ay = (28.6 ± 8.5) x 10 , com- 



pared to an expected slope of (10.0±2.3) x 10 -2 . We use 
the full covariance matrix for the corrected Avb values 
when minimizing x 2 m order to account for the correla- 
tions between bins in the parton-level distribution. The 
systematic uncertainties on Avb in each bin are added 
in quadrature to the diagonals of the covariance matrix. 
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TABLE VIII: The asymmetry at the parton level as measured in the data, compared to the SM tt expectation, as a function 
of |Ay|. 



Parton level 


Data 


SM tt 


|Ay| 


Afb ± stat ± syst 


A FB 


0.0 - 0.5 


0.048 ± 0.034 ± 0.022 


0.023 ± 0.007 


0.5 - 1.0 


0.180 ± 0.057 ± 0.038 


0.072 ± 0.022 


1.0 - 1.5 


0.356 ± 0.080 ± 0.036 


0.119 ± 0.036 


> 1.5 


0.477 ± 0.132 ± 0.070 


0.185 ± 0.056 


< 1.0 


0.101 ± 0.040 ± 0.024 


0.043 ± 0.013 


> 1.0 


0.392 ± 0.093 ± 0.042 


0.139 ± 0.042 



VII. DEPENDENCE OF THE ASYMMETRY 

ON M ti 

The dependence of Ape on the invariant mass of the 
tt system was also studied in the 5 fb -1 analyses [5J f3] 
with only two bins. M t i is correlated with the rapid- 
ity difference Ay, but because Ay depends on the top- 
quark production angle in addition to Af t f, a measure- 
ment of the M t i dependence can provide additional in- 
formation about the underlying asymmetry relative to 
the j4fb(|A?/|) measurement. In the previous publica- 
tions [51 the CDF and DO measurements of Afb 
at small and large M t i were consistent within statisti- 
cal uncertainties but had quite different central values, 
leading to an ambiguity in the comparison of the results 
and their interpretation. We use the full CDF data set 
and the new techniques introduced in this analysis to 
clarify the dependence of Afb on M t i. 



m 
< 



0.6 1 

— CDF Data, 9.4 fb" 1 

q_5 _ = (8.9 ± 2.3) x lO^GeV/c 2 )" 1 

Total Prediction 

0.4 - a K = (2.4 ± 0.6) x 10 4 (GeV/c 2 )- 1 

0.3- 

0.2- 

0.1 - 




350 400 450 500 550 600 650 700 750 800 

M tf (GeV/c 2 ) 



FIG. 16: The reconstruction-level forward-backward asym- 
metry as a function of M t i with a best-fit line superimposed. 
The last bin contains overflow events. The errors on the data 
are statistical, and the shaded region represents the uncer- 
tainty on the slope of the prediction. 



.- 500 
CO 



«2 400 

c 

> 

UJ 300 



200 



100 



CDF Data - Bkg, 9.4 fb" 1 , Ay > 
CDF Data - Bkg, 9.4 fb" 1 , Ay < 



350 400 450 500 550 600 650 700 750 800 

,2 X 



M tf (GeV/cO 



(a) 




350 400 450 500 550 600 650 700 750 800 

M tf (GeV/c 2 ) 

(b) 



17: (a) 



M t i after background sub tract ion in events 
and negative Aj/ and 1(b) background- 



FIG 

with positive and negative Ay and 
subtracted Apb as a function of M t i with a best-fit line 
superimposed. The last bin contains overflow events. Error 
bars include both statistical and background-related system- 
atic uncertainties. The shaded region in (b) represents the 



theoretical uncertainty on the slope of the prediction. 



We divide the data into several mass bins and de- 
termine the number of events with positive (Np ) and 
negative (Nb) Ay in each bin, from which we calculate 



16 



m 

700 
£ 600 

1 500 

UJ 

400 
300 
200 
100 




-1.5 -1 



— CDF Data - Bkg, 9.4 fb 1 
A FB = 0.030 ± 0.031 

ti Prediction 

A„ = 0.01 8 ±0.012 




-0.5 



0.5 



1.5 2 

Ay 



(a) 



m 250 
© 

£ 200 
c 

■2 150 



100 
50 




-1.5 



■ CDF Data - Bkg, 9.4 fb" 1 
A FB = 0.1 97 ±0.043 

tt Prediction 

A„ = 0.062 ±0.012 



fazh-I- 



-0.5 



0.5 



1 1.5 2 

Ay 



(b) 



FIG. 18: The background-subtracted Ay distributions for events with (a) M t j < 450 GeV/c and (b) M tf > 450 GeV/c 



Error bars include both statistical and background-related systematic uncertainties. 



the asymmetry as a function of M t t according to 



Afb(Mu) 



N F (M t t) - N B (M t t) 
N F (Mu)+N B (Mu)- 



(7) 



The M t j-dependent asymmetry is compared to the NLO 
ti plus background prediction in Fig. [16] and Tabic |IX| 
The M t i spectrum is divided into intervals of 50 GeV/c 2 
below 600 GeV/c 2 and 100 GeV/c 2 intervals above 
600 GeV/c 2 , with the final bin containing overflow 
events. The M t i resolution across this range varies as 
a function of mass, being approximately 50 GeV/c 2 at 
the lowest masses and increasing to near 100 GeV/c 2 
at very high mass. A linear fit of the observed data 
has x 2 /Ndof = 1-0/5 and yields a slope of olm u = 
(8.9 ± 2.3) x 10~ 4 (GeV/c 2 )-\ which is non-zero with 
significance in excess of 3<r. The predicted slope at the 
reconstruction level is (2.4 ±0.6) x 10~ 4 (GeV/c 2 )" 1 . 

After removing the background contribution, 
Fig. 17(a) compares the observed M t i distributions 
in forward and backward events, with an excess of 
forward events in many bins. These distributions 
are converted into asymmetries as a function of 
M t £, as shown in Fig. 17(b)| and summarized in 



Table |XJ The linear fit to the background-subtracted 
asymmetries yields x 2 /^dof = 1.1/5 and a slope of 
(10.9 ± 2.8) x 10" 4 (GeV/c 2 )" 1 , with the predicted 
slope being (3.0 ± 0.7) x 10" 4 (GeV/c 2 )- 1 . 

At the background-subtracted level, we divide the 
data into two regions of M t i (above and below 
450 GeV/c 2 ) for direct comparison to the 5.3 fb -1 
CDF analysis [2]. The Ay distributions at high and 
low mass are shown in Fig. |18| yielding asymmetries of 
0.030 ± 0.031 for M a < 450 GeV/c 2 and 0.197 ± 0.043 
for M t i > 450 GeV/c 2 , where the uncertainties include 
statistical and background-related systematic contribu- 



tions. These are in good agreement with the values 
from the 5.3 fb _1 analysis, which found background- 
subtracted asymmetries of —0.022 ± 0.043 for M t t < 
450 GeV/c 2 and 0.266±0.62 for M tl > 450 GeV/c 2 0. 
To check against potential systematic effects, the be- 
havior of the background-subtracted asymmetry at high 
and low M t i in various subsets of the data is summa- 
rized in Table |XI| The M t i dependence is consistent 
across lepton charge and lepton type. It is consistent 
(within relatively large statistical uncertainties) across 
single- and double-&-tagged events. The asymmetry is 
larger in events with exactly four jets than it is in events 
with at least five jets, an effect that is discussed further 
in Sec. El 

We determine the parton-level mass dependence of 
^4fb by correcting the Ay and M t t distributions si- 
multaneously. The correction procedure operates on 
one-dimensional distributions, so we apply the algo- 
rithm to the observable sign(Ay) A/ t j , a one-dimensional 
mass distribution where M t j is reflected across the y- 
axis for events with Ay < 0. The resulting M t i dis- 
tributions are shown in Fig. 19(a) These distribu- 



tions are combined to determine the differential asym- 
metry as a function of M t j shown in Fig. |19(b)| and 
summarized in Table IXIll The best-fit line to the mea- 
sured data asymmetries at parton level has a slope «m iE 
= (15.2 ± 5.0) x 10~ 4 (GeV/c 2 )-\ compared to the 
powheg prediction of (3.4 ± 1.2) x 10~ 4 (GeV/c 2 ) -1 . 



VIII. DETERMINATION OF THE 
SIGNIFICANCE OF THE DEPENDENCE OF 
THE ASYMMETRY ON |Ay| AND M tl 

The slopes of the linear dependencies of Afb on | Ay| 
and Mu provide a measure of the consistency between 
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TABLE IX: The asymmetry observed in the reconstructed data, compared to the SM tt plus background expectation, as a 
function of M t j. 





Data 


SM tt + Bkg. 


Mu (GeV/c 2 ) 


Afb ± stat 


^FB 


< 400 


-0.005 ± 0.030 


0.002 ± 0.006 


400 - 450 


0.053 ± 0.039 


0.017 ± 0.010 


450 - 500 


0.118 ± 0.050 


0.028 ± 0.012 


500 - 550 


0.152 ± 0.067 


0.040 ± 0.018 


550 - 600 


0.128 ± 0.086 


0.067 ± 0.025 


600 - 700 


0.275 ± 0.101 


0.054 ± 0.024 


> 700 


0.294 ± 0.134 


0.101 ± 0.042 



TABLE X: The asymmetry at the background-subtracted level as measured in the data, compared to the SM tt expectation, 
as a function of Mu . 





Data 


SM tt 


Mu (GeV/c 2 ) 


^4fb ± (stat+syst) 


A FB 


< 400 


0.003 ± 0.038 


0.012 ± 0.006 


400 - 450 


0.076 ± 0.049 


0.031 ± 0.011 


450 - 500 


0.149 ± 0.061 


0.039 ± 0.015 


500 - 550 


0.198 ± 0.083 


0.060 ± 0.022 


550 - 600 


0.156 ± 0.104 


0.083 ± 0.030 


600 - 700 


0.361 ± 0.128 


0.077 ± 0.028 


> 700 


0.369 ± 0.159 


0.137 ± 0.049 



TABLE XI: Various measured asymmetries after background subtraction, inclusively and at small and large M t j. 



^4fb ± (stat+syst) 



Sample 


Inclusive 


Mu < 450 GeV/c 2 


M t i > 450 


GeV/ 


All data 


0.087 


± 


0.026 


0.030 


± 


0.031 


0.197 


± 


0.043 


Positive leptons 


0.094 


± 


0.036 


0.034 


± 


0.044 


0.207 


± 


0.060 


Negative leptons 


0.080 


± 


0.035 


0.027 


± 


0.043 


0.186 


± 


0.057 


Exactly 1 b tags 


0.100 


± 


0.031 


0.047 


± 


0.036 


0.220 


± 


0.049 


At least 2 6 tags 


0.037 


± 


0.045 


-0.018 


± 


0.055 


0.134 


± 


0.073 


Electrons 


0.079 


± 


0.039 


0.017 


± 


0.047 


0.195 


± 


0.062 


Muons 


0.094 


± 


0.033 


0.041 


± 


0.040 


0.197 


± 


0.055 


Exactly 4 jets 


0.110 


± 


0.031 


0.029 


± 


0.037 


0.256 


± 


0.049 


At least 5 jets 


0.033 


± 


0.044 


0.034 


± 


0.053 


0.033 


± 


0.077 



the data and the SM prediction. We quantify this con- 
sistency in a more rigorous manner by repeating the 
measurement on large ensembles of simulated experi- 
ments generated according to the SM prediction and 
determining the probabilities, or p-values, for observ- 
ing the actual data given the SM assumption. Each 
p-value is defined as the fraction of simulated experi- 
ments in which the measured slopes are at least as large 
as those found in the data, 



.simulated 



> a 



data 



We use the background-subtracted sample for mea- 
suring these p-values because it provides access to an 
asymmetry calculation that has been corrected for back- 



ground dilution but is still independent of the assump- 
tions needed to extract parton-level information. We 
start from the predicted distribution at the reconstruc- 
tion level, created from the standard model predictions 
of POWHEG and the various background contributions 
proportioned as in Table [TTJ The population of each bin 
of this predicted distribution is fluctuated within its un- 
certainty, which includes the statistical uncertainty on 
the contents of that bin, the systematic uncertainties on 
the various background contributions, as described in 



Sec. VB above, and the theoretical uncertainty on the 
powheg prediction. For each simulated experiment, 
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TABLE XII: The asymmetry at the parton level as measured in the data, compared to the SM tt expectation, as a function 
of M t - t . 



Parton level 


Data 


SM tt 


M tl (GeV/c 2 ) 


Ape ± stat ± syst 


A FB 


< 450 


0.084 ± 0.046 ± 0.026 


0.047 ± 0.014 


450 - 550 


0.255 ± 0.062 ± 0.028 


0.090 ± 0.027 


550 - 650 


0.370 ± 0.084 ± 0.084 


0.117 ± 0.035 


> 650 


0.493 ± 0.158 ± 0.077 


0.143 ± 0.043 


< 450 


0.084 ± 0.046 ± 0.026 


0.047 ± 0.014 


> 450 


0.295 ± 0.058 ± 0.031 


0.100 ± 0.030 
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(a) The parton-level M t i di strib utions for events 
and negative Ay and (b) the parton-level 



positive and negative Ay and 
forward-backward asymmetry as a function of M t i with a 
best-fit line superimposed. The last bin contains overflow 
events. Uncertainties are correlated and include both sta- 
tistical and systematic contributions. The shaded region in 
(b) represents the theoretical uncertainty on the slope of the 
prediction. 



the nominal background prediction with the normaliza- 
tions of Table [TT] is subtracted, and the slopes of the 
remaining asymmetries as a function of |Ay| and M t i 
are fit. We find p-values of 14.7 x 1CT 3 for A FB {\Ay\) 
and 7.4 x 10~ 3 for A FB (M t t), corresponding to 2.2a and 
2.4ct discrepancies respectively (based on a one-sided in- 
tegration of the normal probability distribution). 

We also consider the intercepts associated with the 
linear fits. As shown in the fits of Sections |VI| and [VI I| 
A FB (\Ay\) and A FB (M t i) vanish at |Ay| = and M ti = 
2mtop within good approximation, a behavior that is 
expected. In the case of A F b(|A?/|) in particular, the 
requirement that ^4fb = at \Ay\ = follows directly 
from the continuity and differentiability of the Ay dis- 
tribution. This a priori belief in the smoothness of Ay 
can be used to impose an additional constraint on the fit 
of the observed Afb(|A?/|) distribution. With this con- 
straint that Afb(|A7/| = 0) = 0, the observed slope after 
background subtraction changes from (18.2±5.7) x 10~ 2 
in Fig. pi to (15.5 ± 3.3) x 1CT 2 , with the slope pre- 



dicted by powheg changing from (6.6 ± 1.6) x 10~ 2 to 
(5.3 ±1.0) x 10~ 2 . The p- value for the prediction to 
fluctuate such that it produces a slope at least as large 
as the data decreases from 14.7 x 10~ 3 to 2.2 x 10~ 3 , 
corresponding to an increase in significance from 2.2a 
to 2.8cr. 



IX. DEPENDENCE OF THE ASYMMETRY 
ON THE TRANSVERSE MOMENTUM OF THE 
tt SYSTEM 

The QCD asymmetry at NLO arises from the sum of 
two different effects IT]. The interference of the 2^2 



LO tree-level diagrams (upper left of Fig. 20 1 and the 
NLO box diagrams (upper right) produces a positive 
asymmetry ("Born-box" interference), while the inter- 
ference of 2 — > 3 tree-level diagrams with initial-state 
(lower left) and final-state radiation (lower right) pro- 
duces a negative asymmetry ( "ISR-FSR" interference) . 
In the latter final state, tt plus an additional jet, the 
tt system acquires a transverse momentum , while in 
the former case with an exclusive final state, all events 
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FIG. 20: Interfering qq 
diagrams. 



tt (top) and qq — > itj (bottom) 



have Pj. — 0. The resultant SM asymmetry at NLO is 
therefore the sum of two effects of different sign, with 
very different p" dependence. The virtual effects from 
Born-box interference are larger, leading to a net posi- 
tive asymmetry. Recent work has also emphasized that 
color coherence during the hadronization process can 
produce a significant p" dependence for the asymmetry 
in Monte Carlo generators that include hadronization, 
with the degree of the p" dependence varying greatly 
depending on the details of the implementation of color 
coherence [10 . The verification of the p" dependence 
of the asymmetry is therefore crucial to understanding 
the reliability of the SM predictions for Afb [I] , as well 
as testing for possible new effects beyond the SM. 

In this section, we first compare and discuss several 
predictions for Ap^(p^). We then compare the data to 
two of these predictions (the NLO with hadronization 
prediction from powheg and the LO with hadroniza- 
tion prediction from PYTHIA), showing that the asym- 
metry in the data displays the same trend with respect 
to p% as observed in both POWHEG and pythia, and 
that the excess inclusive asymmetry in the data is con- 
sistent with a p"-independent component. 

We define the p" dependence of the asymmetry as 



0$) 



N F (pf) 



N F (p%) + N B (p%) 



(8) 



The expected SM parton-level asymmetry is shown 
for four predictions in Fig. [21] The matrix ele- 
ments for pythia are LO for tt production, with some 
higher-order effects approximated through hadroniza- 
tion. There is essentially no net inclusive asymmetry 
in pythia due to the underlying 2 — > 2 matrix ele- 
ments in the hard-scattering process; gluon emission 
during hadronization results in a negative asymmetry 
for non-zero p" events, leaving a positive asymmetry in 
the low-p" region. 

The other three curves suggest a different behavior for 
the pj, dependence at NLO. The MCFM calculation uses 
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FIG. 21: Expected Afb as a function of the p" of the tt sys- 
tem at the parton level from mcfm, powheg, and pythia, 
as well as a NLO prediction for events where the top-quark 
pair is produced in association with an extra energetic jet. 



NLO matrix elements for tt production, and includes 
both the Born-box and ISR-FSR interference terms, 
with the result being a parton-level output with two 
partons (tt) or three partons (ft plus a gluon) in the final 
state. In MCFM, events produced by the virtual matrix 
elements with Born-box interference have p" = and a 
positive asymmetry, while events produced by the real 
matrix elements describing gluon radiation have non- 
zero p" and a negative asymmetry, powheg has the 
same NLO matrix elements as mcfm, with additional 
higher-order effects approximated through hadroniza- 
tion performed by pythia. The additional radiation 
from the hadronization process results in a migration 
of events in p" and thus a moderation of the otherwise 
bimodal p" behavior observed in MCFM. 

The POWHEG prediction with pythia hadronization 
can be partially checked against a recent NLO calcula- 
tion for tt production in association with an extra ener- 
geticjet (p£* > 20 GeV/c and |r? je t| < 2.0) [33, shown 
as "tt+jet" . This calculation has a Born-level final state 
with three partons (tt plus a gluon), and thus it is most 
relevant for comparison to the other predictions at high 
p". It contains virtual matrix elements for the tt+jet 
final state as well as real corrections from final states 
with tt and two extra jets. The negative asymmetry ob- 
served in the tree-level prediction for tt+jet (as shown 
in MCFM at high-p") is reduced with the full NLO cal- 
culation of this final state. In the high-p" region, we 
see that the powheg predictions are in good agreement 
with those from the NLO tt+jet calculation. 

In Fig. [3] we show that the reconstructed p% spec- 
trum in the data is well-reproduced by the tt signal 
and background model simulations. Building on this, 
we study the p" d ependence of the asymmetry in the 
data. Figure 22 shows Afb (Pt) f° r the data after 
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FIG. 22: The background-subtracted forward-backward 
asymmetry in the data as a function of the transverse mo- 
mentum of the ti system, compared to both POWHEG and 
pythia. Error bars include both statistical and background- 
related systematic uncertainties. The last bin contains over- 
flow events. 



FIG. 24: The background-subtracted forward-backward 
asymmetry in the data as a function of the transverse mo- 
mentum of the ti system, compared to both POWHEG and 
pythia. The model predictions have been normalized by 
the addition of AAnlo to powheg and AAlo to pythia 
as described in the text. Error bars include both statistical 
and background-related systematic uncertainties. The last 
bin contains overflow events. 
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FIG. 23: The forward-backward asymmetry as a function of 
the transverse momentum of the ti system for three models 
at the background-subtracted level: Octet A, SM pythia, 
and SM pythia normalized by the addition of AAo ct . The 
last bin contains overflow events. 



background subtraction compared to predictions from 
POWHEG (hadronized with pythia) and from PYTHIA. 
The bin size is equal to one-half of the detector reso- 
lution in p". The trends of the parton- level curves in 
Fig. [21] are reproduced: the LO prediction has a steady 
drop, while the NLO prediction tends to zero or slightly 
below. The data show a clear decrease with p", but lie 
above the models. We investigate this using the ansatz 
that the data contain an additional source of asymme- 
try that is independent of p". In this case, because 
independent asymmetries are additive, it should be pos- 
sible to normalize the model predictions to the data by 
adding a constant offset AA that is equal to the excess 
observed inclusive asymmetry in the data. 



We test this ansatz using the color-octet model 
Octet A (implemented in madgraph with hadroniza- 
tion performed by pythia) described at the end of 
Sec.JTlj In this LO model, the octet physics induces 
a ^-independent inclusive ti asymmetry 0.106 at the 
background-subtracted level (we neglect very small sta- 
tistical uncertainties in these large Monte Carlo sam- 
ples). We wish to compare the dependence of 
this asymmetry to the LO PYTHIA model, which has 
a background-subtracted asymmetry of —0.021. The 
inclusive difference is A^4oct = 0.127. If the excess 
asymmetry in Octet A is independent of p", we ex- 
pect that A^™ a (pt) + AA 0ct reproduces satisfactorily 
Figure ! 



A Octet A-t^tt 
^FB \Pt 



2.3 



shows this test in the simulated 
samples, with the p"-dependent behavior of Octet A 
being described well by the addition of the constant 
normalization factor AA 0c t to A^g HIA (p"). 

We use this procedure to normalize the ^4fb(p") 
models of powheg and pythia to the total inclusive 
asymmetry observed in the data. Since this artificial 
procedure adjusts the mean values such that they are 
exactly equal, we do not assign uncertainties to the off- 
sets. The asymmetry after background subtraction is 
0.087 in the data, 0.033 in NLO powheg (Table [III]), 
and —0.021 in LO pythia, resulting in offset terms 
AA NLO = 0.054 and AA LO = 0.108. 

The normalized ^4fb(Pt) m odels are compared to the 
data in Fig. [24[ Within the experimental uncertainties, 
the A FB (p") behavior of the data is described well by 
both models. We conclude that the reconstruction and 
modeling of the dependence of the asymmetry is 
robust, and that the excess asymmetry in the data is 
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FIG. 25: The measured asymmetries as a function of (a) \ Ay\ and (b) M t j in the data at the three different levels of correction. 
The error bars include both statistical uncertainties and the appropriate systematic uncertainties for each correction level 
as described in the text. The last bin contains overflow events. 



TABLE XIII: The measured inclusive forward-backward asymmetry and the best-fit slopes for Afb(|Ai/|) and ApB(A^tt) 
at the different levels of correction. The uncertainties include the statistical uncertainties and the appropriate systematic 
uncertainties for each correction level as discussed in the text. 



Correction level 


Inclusive 
A FB 


Slope 


Slope 

aM ti 


Reconstruction 
Background-subtracted 
Parton 


0.063 ± 0.019 
0.087 ± 0.026 
0.164 ± 0.045 


(14.3 ±4.3) x 10~ 2 
(18.2 ±5.7) x 10" 2 
(28.6 ±8.5) x 10" 2 


(8.9 ±2.3) x 10" 4 (GeV/c 2 )- 1 
(10.9 ±2.8) x KT 4 (GeV/c 2 )" 1 
(15.2 ±5.0) x KT 4 (GeV/c 2 )- 1 



consistent with being independent of p". 

Finally we note the connection between and jet 
multiplicity. In events with one or more extra ener- 
getic jets, we expect the ti system to have large 
due to recoil against these additional jets. In Table |Xl| 
and Ref. [5] a difference was noted in the asymmetry 
measurements at the background-subtracted level be- 
tween events with exactly four jets_ and at least five 
jets. Rephrasing this in terms of p", we find that the 
mean in five-jet events is 34.4±0.6 GeV/c compared 
to 18.6 ± 0.3 GeV/c in events with only four jets. The 
smaller asymmetry in events with extra jets is seen to 
be consistent with the observed A F b(Pt) behavior. 

X. CONCLUSIONS 

We study the forward-backward asymmetry A F b in 
top-quark pair production using the full CDF Run II 
data set, finding results that are consistent with previ- 
ous measurements with reduced uncertainties. In the 
reconstruction- level data, prior to any corrections, we 
observe an asymmetry of 0.063 ± 0.019(stat) consistent 
with a linear dependence on both \Ay\ and M t i. The 



inclusive asymmetry exceeds by approximately 2a the 
prediction of 0.020 ± 0.012 that results from the sum of 
the NLO standard model as implemented in powheg 
and the background expectation. 

After removing the predicted background contribu- 
tion, the inclusive asymmetry is larger and the approxi- 
mately linear behavior of A F b with increasing |Ay| and 
M t i remains. We determine the significance of the ra- 
pidity and mass dependence by comparing the best-fit 
slopes in the data to the POWHEG prediction, finding p- 
values of 14.7 x 10~ 3 for A FB as a function of |Ay| and 
7.4 x 10 -3 for Afb as a function of M t j, corresponding 
to significances of 2.2c and 2.4er respectively. 

Finally, we correct the results to the parton level to 
find the differential cross section da/d{Ay) and directly 
compare with theoretical predictions, finding an inclu- 
sive parton-level asymmetry of 0.164±0.045(stat±syst). 
We measure a linear mass dependence A F b (M t j) at par- 
ton level with slope (15.2 ± 5.0) x 10~ 4 (GeV/c 2 )- 1 , 
exceeding the NLO standard model prediction of (3.4 ± 
1.2) x 10~ 4 (GeV/c 2 )" 1 by 2.3cr, and a linear rapidity 
dependence A FB (\Ay\) with slope (28.6 ± 8.5) x 10~ 2 , 
exceeding the prediction of (10.0 ± 2.3) x 10~ 2 by 2.1cr. 

The measured functional dependencies of the asym- 
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metry at the various stages of the correction procedure 
are compared directly in Fig. [25] The inclusive A FB 
values at each level and the measured slopes of the \ Ay\ 
and M t i dependencies of the asymmetry are compared 
in Table Knil 

In addition to these measurements, we also find that 
the dependence of Afb on the transverse momentum of 
the tt system is in good agreement with the SM pre- 
diction, under the assumption that the excess observed 
asymmetry is independent of . 

The larger-than-expected top-quark forward- 
backward asymmetry observed in previous analyses 
at the Tevatron persists in the full CDF II data set. 
The results reported here are consistent with previous 
measurements of the asymmetry. A detailed study 
shows linear dependencies of the asymmetry as a func- 
tion of the invariant mass of the top-quark pair and 
the rapidity difference |Ay| that are similar to what 
is expected from the SM but with larger slopes than 
current calculations predict. This new measurement 
serves as a means to understand better the possible 
higher-order corrections to the SM prediction as well 
as potential effects from non-SM physical processes. 
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Appendix: Validation of charge asymmetry 
measurements with the CDF II detector 



With a pp initial state and a tt final state that 
are symmetric under charge conjugation, the forward- 
backward asymmetry is equivalent to a charge asym- 
metry. The asymmetry measurements rely crucially 
on measurement of the lepton charge to determine the 
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FIG. 26: The lepton pr distribution in events with a W bo- 
son and only one observed jet. The last bin contains overflow 
events. 



charges of all reconstructed particles in the tt final state. 
This is particularly important when the lepton px is 
large, as large M t i is correlated with large-lepton-p-r 
events, in which the determination of the lepton charge 
is more challenging. It is therefore important to verify 
that the lepton charge is measured correctly over the 
lepton pt and rj ranges pertinent to the tt measure- 
ments. 

We do this in the large sample of CDF events con- 
taining a W boson and only one observed hadronic 
jet. In addition to an abundant, low-background sig- 
nal, this sample features a well-understood, lepton-p-r- 
dependent asymmetry in the direction of motion of the 
lepton from the tU-boson decay, which is used to gauge 
the charge measurement. We measure the asymmetry 
in the observable q-r]i ep , where q is the lepton charge and 
7]i op is the pseudorapidity of the lepton. At low lepton 
Pt, the asymmetry is positive and dominated by asym- 
metric contributions to the proton parton-distribution 
function from u and d quarks, while at large lepton 
Pt, the asymmetry is negative and dominated by ef- 
fects from the electroweak decay of the W boson. We 
compare the data and prediction for the leptonic asym- 
metry over the relevant ranges of lepton px and r7i ep to 
test whether we reproduce the known SM asymmetries 
in this important control region. As in the analysis of 
the tt signal sample, SM JU+jet production is modeled 
using the ALPGEN generator [27] . 

The lepton selection in this sample is the same as that 
for tt candidate events. We require there to be only one 
observed hadronic jet, for which no 6-tag requirement 
is applied. We also release the Ht requirement for this 
sample. Finally, we introduce a new variable, the mini- 
mum VF-boson mass. We add the four-momenta of the 
identified lepton and the "neutrino" , which is defined to 
be a massless particle with the x- and y-components of 
momentum given by the $t, and the z-component cho- 



23 




FIG. 27: The forward-backward asymmetry in q ■ r]\ cp at the detector level as a function of (a) |?7i C p| and (b) lepton px in 
events with a W boson and only one observed jet. The error bars on the data include only a statistical contribution, with 
the uncertainty on the SM prediction shown as a band around the predicted asymmetry. The last bin contains overflow 
events. 



TABLE XIV: The q ■ r/i C p asymmetry in the W+l jet sample, compared to SM expectations, for small and large lepton px- 





Lepton pr 


< 60 GeV/c 


Lepton pr > 60 GeV/c 


Observed data 


0.083 


± 0.001 


-0.009 ± 0.004 


SM prediction 


0.089 


± 0.004 


-0.001 ± 0.013 


Data minus prediction 


-0.006 


± 0.004 


-0.008 ± 0.014 



TABLE XV: The q ■ r/i cp asymmetry in the W+l jet sample, compared to SM expectations, for small and large |»7i ep |. 







< 0.75 


k/icp > 0.75 


Observed data 


0.059 


± 


0.001 


0.124 


± 


0.002 


SM prediction 


0.063 


± 


0.005 


0.134 


± 


0.008 


Data minus prediction 


-0.004 


± 


0.005 


-0.010 


± 


0.008 



sen to minimize the total mass of the lepton+neutrino 
system. We require this mass to exceed 20 GeV/c 2 , 
removing most of the non-VF contribution to this data 
sample. After applying this selection, we have approxi- 
mately 800 000 total data events. 

The lepton px in the W+l jet sample is given in 
Fig. [26j Good agreement with the prediction is seen 
over the entire spectrum. Compared to the lepton px in 
tt decays shown in Fig.[5j this distribution is softer, but 
it still provides sufficient precision in the high-lepton-pT 
region relevant to the it sample. Figure 27 shows the 
asymmetries in q ■ ri\ cp as a function of |ryi ep | and lepton 
Pt- Across the entire spectrum, the asymmetry is mea- 
sured with good accuracy and is in excellent agreement 
with the SM prediction. The biggest difficulty with the 



comparison is the uncertainty due to model sampling at 
very large |j?iep| and lepton px- 

In ti events, lepton px is correlated with M t j, with 
higher mass tt pairs leading to larger lepton px ■ There- 
fore, in the context of Ayb {M t +), where a large asym- 
metry is observed at high mass, we are particularly in- 
terested in events with high lepton px- The measured 
asymmetries in two bins of lepton px are given in Ta- 



ble XIV for direct comparison to the SM prediction. In 
the context of Afb(|A?/|), we also list the asymmetries 
for two bins of | r7i ep | in Table XV Excellent agreement, 
to within 1%, is found between the data and the pre- 
diction in all regions of lepton px and ryiep using this 
high-precision control sample, supporting confidence in 
the determination of the lepton charge. 
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